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ABSTRACT 
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making the Hartree approximation, although the numerical coefficients differ. 



PACS numbers: 04.30.Nk, 04.62. +v, 98.80. Cq, 98.80.Hw 

' e-mail: miao@phys.ufl.edu 
* e-mail: woodard@phys.ufl.edu 



1 Introduction 



Gravitons and massless, minimally coupled scalars can mediate vastly en- 
hanced quantum effects during inflation because they are simultaneously 
massless and not conformally invariant pp. One naturally wonders how in- 
teractions with these quanta affect themselves and other particles. The first 
step in answering this question on the linearized level is to compute the one 
particle irreducible (1PI) 2-point function for the field whose behavior is in 
question. This has been done at one loop order for gravitons in pure quantum 
gravity [2j , for photons [3], H] and charged scalars [5] in scalar quantum elec- 
trodynamics (SQED), for fermions [Bl El and Yukawa scalars [8] in Yukawa 
theory, for fermions in Dirac + Einstein [9] and, at two loop order, for scalars 
in 4 theory [TDj. The next step is using the 1PI 2-point function to correct 
the linearized equation of motion for the field in question. That is what we 
shall do here for the fermions of massless Dirac + Einstein. 

It is worth reviewing the conventions used in computing the fermion self- 
energy [pj. We worked on de Sitter background in conformal coordinates, 

ds 2 = a 2 (r])(—dri 2 + dx-dx] where a(n) = = e Ht . (1) 

v ' Hi] 

We used dimensional regularization and obtained the self-energy for the con- 
formally re-scaled fermion field, 

= a ( VV( x ) • (2) 

The local Lorentz gauge was fixed to allow an algebraic expression for the 
vierbein in terms of the metric [UJ . The general coordinate gauge was fixed 
to make the tensor structure of the graviton propagator decouple from its 
spacetime dependence [12j [13] . The result we obtained is, 
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where k 2 = lG^G is the loop counting parameter of quantum gravity. The 
various differential and spinor-differential operators are, 

d 2 = rrdpd v , V 2 ee dA , ? ee 7 ^ and ^ ee 7 % , (4) 
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where if is the Lorentz metric and 7^ are the gamma matrices. The con- 
formal coordinate interval is basically Ax 2 = (x — — x') u i] t j, u , up to a 
subtlety about the imaginary part which will be explained shortly 
The linearized, effective Dirac equation we will solve is, 



In judging the validity of this exercise it is important to answer five questions: 

1. What is the relation between the (D-number, effective field equation (GO) 
and the Heisenberg operator equations of Dirac + Einstein? 

2. How do solutions to (J5]) change when different gauges are used? 

3. How do solutions to (JSJ) depend upon the finite parts of counterterms? 

4. What is the imaginary part of Arc 2 ? and 

5. What can we do without the higher loop contributions to the fermion 
self-energy? 

Issues 1 and 2 are closely related, and require a lengthy digression that we 
have consigned to section 2 of this paper. In this Introduction we will com- 
ment on issues 3-5. 

Dirac + Einstein is not perturbatively renormalizable [14] . so we could 
only obtain a finite result by absorbing divergences in the BPHZ sense [15} 
[16j [TTJ, [18] using three higher derivative counterterms, 



No physical principle seems to fix the finite parts of these counterterms so 
any result which derives from their values is arbitrary. We chose to null local 
terms at the beginning of inflation (a = 1), but any other choice could have 
been made and would have affected the solution to (j3J). Hence there is no 
point in solving the equation exactly. However, each of the three countert- 
erms is related to a term in ([3]) which carries a factor of In (aa'), 





2 H 2 \-^—$d 2 + a 2 D(D-l)$+a 3 $}8 D (x-x') . 



(6) 




(8) 
(9) 




2 



Unlike the aj's, the numerical coefficients of the right hand terms are uniquely 
fixed and completely independent of renormalization. The factors of ln(aa') 
on these right hand terms mean that they dominate over any finite change in 
the aj's at late times. It is in this late time regime that we can make reliable 
predictions about the effect of quantum gravitational corrections. 

The analysis we have just made is a standard feature of low energy ef- 
fective field theory, and has many distinguished antecedents P3H [201 ED E2J 
E3l [2H [251 [2HI I2ZI [211 [2HI SOI EH [321 S3]. Loops of massless particles make 
finite, nonanalytic contributions which cannot be changed by counterterms 
and which dominate the far infrared. Further, these effects must occur as 
well, with precisely the same numerical values, in whatever fundamental the- 
ory ultimately resolves the ultraviolet problems of quantum gravity. 

We must also clarify what is meant by the conformal coordinate interval 
Ax 2 (x;x') which appears in ([3]). The in-out effective field equations corre- 
spond to the replacement, 

Ax 2 (x; x 1 ) — ► Ax 2 ++ (x; x') = \\x — x \\ 2 — (| r\ — r( \ —i5) 2 . (10) 

These equations govern the evolution of quantum fields under the assumption 
that the universe begins in free vacuum at asymptotically early times and 
ends up the same way at asymptotically late times. This is valid for scattering 
in flat space but not for cosmological settings in which particle production 
prevents the in vacuum from evolving to the out vacuum. Persisting with 
the in-out effective field equations would result in quantum correction terms 
which are dominated by events from the infinite future! This is the correct 
answer to the question being asked, which is, "what must the field be in order 
to make the universe to evolve from in vacuum to out vacuum?" However, 
that question is not very relevant to any observation we can make. 

A more realistic question is, "what happens when the universe is released 
from a prepared state at some finite time and allowed to evolve as it will?" 
This sort of question can be answered using the Schwinger-Keldysh formalism 
[Ml ESI GEE EZ1 ESI ESI HD]- For a recent derivation in the position-space 
formalism we are using, see [JT]. We confine ourselves here to noting four 
simple rules: 

• The endpoints of lines in the Schwinger-Keldysh formalism carry a ± 
polarity, so every n-point 1PI function of the in-out formalism gives 
rise to 2™ 1PI functions in the Schwinger-Keldysh formalism; 



3 



The linearized effective Dirac equation of the Schwinger-Keldysh for- 
malism takes the form (J5J) with the replacement, 



(x; x ) -+- iYlj (x; x ) 5 



(11) 



The ++ fermion self-energy is ((3]) with the replacement (fTUj) ; and 
The +- fermion self-energy is, 
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+ 0(k 4 ),(12) 



with the replacement, 

Ax 2 (x; x') — > Ax 2 _(x; x') = ||x — x^ 2 — (rj — rf + z5) 2 . (13) 

The difference of the ++ and +- terms leads to zero contribution in (JSJ) unless 
the point x' M lies on or within the past light-cone of x M . 

We can only solve for the one loop corrections to the field because we lack 
the higher loop contributions to the self-energy. The general perturbative 
expansion takes the form, 

00 00 
V(x) = K 2e ^ e (x) and [E] (x; x') = £ k u [E e ] (x; x') . (14) 
1=0 e=i 

One substitutes these expansions into the effective Dirac equation fl5]) and 
then segregates powers of k 2 , 



i$Mf°(x) = 



[x:x')ty°(x" 



et cetera. 



^ (15) 

We shall work out the late time limit of the one loop correction *&}(t], x; k, s) 
for a spatial plane wave of helicity s, 
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where k e ^Uj(k, s) = k^Uj(k, s) . 

(16) 
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In the next section we derive the effective field equation. In section 3 we 
derive some key simplifications. In section 4 we solve for the late time limit of 
^fj(r], x; k, s). The result takes the surprising form of a time dependent field 
strength renormalization of the tree order solution. In section 5 we show that 
this can be understood qualitatively using mean field theory. Our results are 
summarized and discussed in section 6. 

2 The Effective Field Equations 

The purpose of this section is to elucidate the relation between the Heisenberg 
operators of Dirac + Einstein — ipi(x), i>i{%) and h^ u (x) - - and the (D- 
number plane wave mode solutions ^i(x; k, s) of the linearized, effective Dirac 
equation (jSJ). After explaining the relation we work out an example, at one 
loop order, in a simple scalar analogue model. Finally, we return to Dirac + 
Einstein to explain how *&i(x; k, s) changes with variations of the gauge. 

2.1 Heisenberg operators and effective field equations 

The invariant Lagrangian of Dirac + Einstein in D spacetime dimensions is, 

C = i^(^-p- 1 )P- 2 )^ 2 )V^+#V(^-^ C ^ cd )V'^ • (17) 

Here e M 6 is the vierbein field and g^ v = e^vcf] hc is the metric. The metric 
and vierbein-compatible connections are, 

(18) 

The Ricci scalar is, 

R = sT (r %,p - r V + r V r V - r Vr CT PM ) . (19) 

The gamma matrices 7^ have spinor indices i,j G {1,2,3,4} and obey the 
usual anti-commutation relations, 

{ 7 b ,j c } = -2 V bc I. (20) 

The Lorentz generators of the bispinor representation are, 

J hc = - A \l\l c }. (21) 
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We employ the Lorentz symmetric gauge, = e^, which permits one 
to perturbatively determine the vierbein in terms of the metric and their 
respective backgrounds (denoted with overlines) [TT] . 

e M = (vW 1 ) " Cub • (22) 
We define the graviton field in de Sitter conformal coordinates as follows, 

9fiu(x) = a 2 (r/> + Kh^xfj where a = . (23) 

By convention the indices of h^ u are raised and lowered with the Lorentz 
metric. We fix the general coordinate freedom by adding the gauge fixing 
term, 

L G¥ = - 1 -a D - 2 r l ^F il F v where F, = tT( Kp,°~\ h p ^ + {D-2)Hah, p 5l) . 

(24) 

One solves the gauge-fixed Heisenberg operator equations perturbatively, 

M*0 = h%(x) + Khl„(x) + K 2 hl u (x) + ... , (25) 
ipi(x) = 4>°(x) + ^1(x) + k 2 iIj 2 (x) + ... . (26) 

Because our state is released in free vacuum at t = (77 = —1/H), it makes 
sense to express the operator as a functional of the creation and annihilation 
operators of this free state. So our initial conditions are that and its 
first time derivative coincide with those of h°^ v (x) at t = 0, and also that 
ipi(x) coincides with ipi(x). The zeroth order solutions to the Heisenberg 
field equations take the form, 

h%(x) = /|pE{v(^;M)e ife «(M) 

+e; u ( m k,\)e- i ^(k,\)}, (27) 

+-^= Vi (k,X)e- ik ^(k,s)} . (28) 

The graviton mode functions are proportional to Hankel functions whose 
precise specification we do not require. The Dirac mode functions Ui(k, s) and 
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Vi(k, s) are precisely those of flat space by virtue of the conformal invariance 
of massless fermions. The canonically normalized creation and annihilation 
operators obey, 

a(k, A), a<(k', A')] = 5 X y{27r) D - 1 5 D - 1 {k-k') , (29) 
{b(ie,8),tf(ft,8')} = 5 ssl {2-a) D - 1 5 D ~\k-k') = [c{k,s),c\k' ,s'j) . (30) 

The zeroth order Fermi field ip®{x) is an anti-commuting operator whereas 
the mode function ^°(x; k, s) is a (D-number. The latter can be obtained from 
the former by anti-commuting with the fermion creation operator, 

—ikr) _ 

^(x;k,s) = a—{^(x),b^k,s)} = Ui (k,s)e^ . (31) 



2k 

The higher order contributions to ipi(x) are no longer linear in the creation 
and annihilation operators, so anti-commuting the full solution ipi(x) with 
b*(k, s) produces an operator. The quantum-corrected fermion mode function 
we obtain by solving (jSJ) is the expectation value of this operator in the 
presence of the state which is free vacuum at t — 0, 

%(x;k,s) = a^faltyiix^tffastykl) . (32) 

This is what the Schwinger-Keldysh field equations give. The more familiar, 
in-out effective field equations obey a similar relation except that one defines 
the free fields to agree with the full ones in the asymptotic past, and one 
takes the in-out matrix element after anti-commuting. 



2.2 A worked-out example 

It is perhaps worth seeing a worked-out example, at one loop order, of the 
relation (|32l) between the Heisenberg operators and the Schwinger-Keldysh 
field equations. To simplify the analysis we will work with a model of two 
scalars in flat space, 

C = -d^Pip - mW - Ax^V-^X^X • (33) 

In this model (p plays the role of our fermion ipi, and x plays the role of 
the graviton h^ u . Note that we have normal-ordered the interaction term 
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to avoid the harmless but time-consuming digression that would be required 
to deal with X developing a nonzero expectation value. We shall also omit 
discussion of counterterms. 

The Heisenberg field equations for (1331) are, 



d 2 X~ \:ip*<p: 
(d 2 - m 2 )ip - \\if 



0, 
. 



As with Dirac + Einstein, we solve these equations perturbatively, 

X {x) = x\x) + \ X \x)+\ 2 x 2 {x) + ... , 
(p(x) = <p°(x) + \ip\x) + X 2 cp 2 (x) + . . . . 

The zeroth order solutions are, 

d D ~ l k re~ ikt 



X\x) 
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(35) 



(36) 
(37) 



(38) 
(39) 



Here k = \\k\\ and u = \Jk 2 + m 2 . The creation and annihilation operators 
are canonically normalized, 

a(k),oJ(k')] = [&(£),&*(#)] = [c{k),c\k')\ = (2n) D - 1 5 D - 1 (k-k') . (40) 

We choose to develop perturbation theory so that all the operators and their 
first time derivatives agree with the zeroth order solutions at t — 0. The first 
few higher order terms are, 
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The commutator of ip°(x) with b*(k) is a C-number, 
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(44) 



However, commuting the full solution with tf(k) leaves operators, 



cp(x), b\k)] = k) + A f dt'f d°- l x' (x 



d 2 — m 2 



x') x (x')<S> (x'-k) 

I ret 
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x\x'),b\k)\ V °(x') + X \x')^{x'-k) 
+X°(x')[^(x'),b\k)}\+0(\ 3 ). (45) 



The commutators in (j4"5l) are easily evaluated 
''dt"! d D ~ l x" lx' ' 



d 2 



X °{x')[p\x'),b\k) 
rt 1 



dt" d B - x x" x' 



x") <p*(x")<p\af)$ u (x";k) , (46) 
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x") x"(x') X "(x")^(x";k). (47) 

/ ret 



Hence the expectation value of (1451) gives, 



(p\ [p(x), b\k)\ \n) = <S>°(x; k) + \ 2 j dt'f cP~V (x 
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x f"dt"f c/ D - 1 x ,, ((x , |^|x ,, ) ict (fi|^ *(x ,, )^ (x , )|fi) 

xVkV) n)WV; £) + 0(A 4 ) . (48) 



To make contact with the effective field equations we must first recognize 
that the retarded Green's functions can be written in terms of expectation 
values of the free fields, 
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(49) 

cm , (so) 

(51) 



i6(t'-t")\ (tt\<p (x')p *(x")\tt) - ({l\<p* {x")<p°{x')\tt) 1 . (52) 



Substituting these relations into 
expression we have been seeking, 



and canceling some terms gives the 
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n 



cp *(x")(p°(x') n) \$°(x"; k) + 0(A 4 ) . (53) 



We turn now to the effective field equations of the Schwinger-Keldysh 
formalism. The (D-number field corresponding to <f(x) at linearized order is 
<&(x). If the state is released at t — then the equation $(a;) obeys is, 

{d 2 - m 2 )$(x) - Jdt'J d D - l x'{M 2 ++ (x\ x) + M 2 _(x; x')}$(x') = . (54) 
The one loop diagram for the self-mass- squared of <p is depicted in Fig. 1. 



x 



x 



Fig. 1: Self-mass-squared for (p at one loop order. Solid lines stands for ip 
propagators while dashed lines represent \ propagators. 

Because the self-mass-squared has two external lines, there are 2 2 = 4 polar- 
ities in the Schwinger-Keldysh formalism. The two we require are [H UT] , 
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h + 0(A 4 ) , (55) 
x 1 ) +0(A 4 ) . (56) 



To recover (|53|) we must express the various Schwinger-Keldysh propaga- 
tors in terms of expectation values of the free fields. The ++ polarity gives 
the usual Feynman propagator |41j . 
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n) . (58) 



The +- polarity propagators are 

i 
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(59) 
(60) 



Substituting these relations into (I55H56I) and making use of the identity 1 
6(t-t?)+6(t'-t) gives, 



Ml + (x; x') + M 2 + _{x] x') = -i\ 2 9(t-t')l (p\ X (x) X ° (x') 



n 



^ (x)^ Q *(x / )|^)-(^|x (^)x (^)|«)(^|/*(^)^(^)|^) +0(A 4 ). (61) 



We now solve (1541) perturbatively. The free plane wave mode function 
(jSj) is of course a solution at order A . With ( |6T1) we easily recognize its 
perturbative development as, 



$(x; k) = $°(x; k) - i\ 2 J^dt'J dP^x' (x 



Q2_ m 2 
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ret 



*V)xV) 



X V)xV)|^)(^WV')|^) 



<p *{x")ip {x') n) }^(x"; k) + 0(A 4 ) . (62) 



That agrees with (I5"3"j) . so we have established the desired connection, 



®(x;k) = (tt \<p{x),tf(k) 



(63) 



at one loop order. 



2.3 The gauge issue 

The preceding discussion has made clear that we are working in a particular 
local Lorentz and general coordinate gauge. We are also doing perturbation 
theory. The function ^®(x; k, s) describes how a free fermion of wave number 
k and helicity s propagates through classical de Sitter background in our 
gauge. What ^}(x;k, s) gives is the first quantum correction to this mode 
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— * 

function. It is natural to wonder how the effective field ^i(x; k, s) changes if 
a different gauge is used. 

The operators of the original, invariant Lagrangian transform as follows 
under diffeomorphisms (x M — > x'^) and local Lorentz rotations (Ay)Q 

ifcx) = K 3 (x'-\x))^(x'-\x)) , (64) 

e'^x) = —Af^xjj^^izj) . (65) 



The invariance of the theory guarantees that the transformation of any solu- 
tion is also a solution. Hence the possibility of performing local transforma- 
tions precludes the existence of a unique initial value solution. This is why no 
Hamiltonian formalism is possible until the gauge has been fixed sufficiently 
to eliminate transformations which leave the initial value surface unaffected. 

Different gauges can be reached using field- dependent gauge transforma- 
tions [12]. This has a relatively simple effect upon the Heisenberg operator 
ipi(x), but a complicated one on the linearized effective field ^i(x; k, s). Be- 
cause local Lorentz and diffeomorphism gauge conditions are typically spec- 
ified in terms of the gravitational fields, we assume x fti and A^- depend upon 
the graviton field h^ u . Hence so too does the transformed field, 

i>'M(x) = A^x'-'Mx^^x'-'Mx)) . (66) 

In the general case that the gauge changes even on the initial value surface, 
the creation and annihilation operators also transform, 

b'[h](k,s) = -^=u*(k,s) f d^'xe-^'M^x) , (67) 
v 2k J 



where rji = — 1/H is the initial conformal time. Hence the linearized effective 
field transforms to, 

%{x;k,s) =a^({l\{tlj' l {h]{x),b'l{h]{k,s)}\n) . (68) 



Of course the spinor and vector representations of the local Lorentz transformation 
are related as usual, with same parameters tu c d(x) contracted into the appropriate repre- 
sentation matrices, 

Aij = Sij ~~ l -u cd J cd l:j + ... and A b c = S b c - cu b c + . . . . 
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This is quite a complicated relation. Note in particular that the h^ v de- 
pendence of x'^[h] and A^f/i] means that ^(x; k, s) is not simply a Lorentz 
transformation of the original function ^i(x;k, s) evaluated at some trans- 
formed point. 



3 Some Key Reductions 

The purpose of this section is to derive three results that are used repeat- 
edly in reducing the nonlocal contributions to the effective field equations. 
We observe that the nonlocal terms of ([3]) contain 1/Ax 2 . We can avoid 
denominators by extracting another derivative, 

-i- = — ln(Ax 2 ) and = d l h n 2 (Ax 2 ) - 2 ln(Ax 2 )l . (69) 

Ax 2 4 Ax 2 8 L J 

The Schwinger-Keldysh field equations involve the difference of ++ and +- 
terms, for example, 

ln(/i 2 Ax 2 + ) ln0 2 Ax 2 J 
Ax 2 Ax 2 




ln 2 (/i 2 Ax 2 + ) - 21n(//Ax 2 + ) - In V Ax 2 J + 2 ln^Ax 2 J . (70) 



We now define the coordinate intervals Arj = rj — rj' and Ax = ||x — x'\\ in 
terms of which the ++ and +- intervals are, 

Ax 2 + = Ax 2 - (|A??|-^) 2 and Ax 2 _ = Ax 2 - (Ar]+i5) 2 . (71) 

When rf > r\ we have Ax 2 + = Ax 2 _ , so the ++ and +- terms in (ITUI) cancel. 
This means there is no contribution from the future. When r]' < r] and 
Ax > Ar) (past spacelike separation) we can take 5 = 0, 

ln(/i 2 Ax 2 + ) = ln[ / u 2 (Ax 2 -Ar/ 2 )] = ln(^ 2 Ax 2 _) (Ax > Ar] > 0) . (72) 

So the ++ and +- terms again cancel. Only for rf < r\ and Ax < Ari (past 
timelike separation) are the two logarithms different, 

ln(/i 2 Ax 2 J = ln[/i 2 (Ar/ 2 -Ax 2 )] ± m [Ar] > Ax > 0) . (73) 
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Hence equation (1701) can be written as 
ln(/i 2 Ax 2 + ) ln(/i 2 Ax 2 _) m 



Ax 2 



Ax 2 



— d 2 \6(Ar]-Ax) \n(/j 2 (Ar] 2 -Ax 2 )-l \.(74) 



This step shows that the Schwinger-Kledysh formalism is causal. 

To integrate (1?4"]) up against the plane wave mode function (|T6|) we first 
pull the x M derivatives outside the integration, then make the change of 
variables x' = x+r and perform the angular integrals, 



(Fx 



, f ln(/i 2 A 
\ Ax\ 
i2tt 2 



xA 



*?{rf,x, k,s) 



Axi 



Ui (k,s)d 2 e ik - X dr/—= drr sin(fcr)| ln[/i 2 (A^ 2 -r 2 )] - 1} 
Jm \2kh L J 



k 
%2n 2 



e ik -* Ui (k,s)[-d 2 -k 2 } I 'drj'e- ikr >' Ar) 



2k Jo 

n 



kV2k 



i e -ikr! , 



x ^Lzsin(az){ln(l-z 2 ) + 21n(^)-l}. (75) 



Here a = kAr] and rji = — 1/H is the initial conformal time, corresponding to 
physical time t = 0. The integral over 2 is facilitated by the special function, 



1 2 1 

dzzsmiaz) ln(l — z 2 ) = — sin (a) - 

a 2 a 2 



cos(a) + asin(a 



x 



si(2a)- 



7T 



+ 



sin (a —a cos (a 



ci(2a)- 7 -ln(-)j . (76) 



Here 7 is the Euler-Mascheroni constant and the sine and cosine integrals 
are, 



si x 



00 , sin(t) 
dt — — 
t 



7T 



h / dt 



sint 







t 



f°° , cos t , , . , , 

ci(x) = - I dt — = 7 + ln(x) + / (it 







cos(t) — 1 



(77) 
(78) 



After substituting the £ function and performing the elementary integrals, 
f5l) becomes, 



d x 



, f ln(/i 2 A 



\ Ax 2 



ln(// 2 Aa^_) 



$°(r/,x,M 



Ax 2 , 



x(9 2 ?? + l)/^'e 



>e,-iW , 



a 2 £(a)+[21n(^)-l 



227T 2 



Ak-x 



kV2k 



Ui(k,s) 



sin a —a cos a 



•■ (79) 
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One can see that the integrand is of order a 3 ln(a) for small a, which 
means we can pass the derivatives through the integral. After some rear- 
rangements, the first key identity emerges, 



d 4 x 



, f hi(/i 2 A 



\ Ax* 



ln(/z 2 Ax+_ 



Ax 2 



f; k, s) 



i4ir 2 k * (77, x; k, s) drj'e 



I ikA-q 



cos(kAri) / dt — — 



+ sin(kAr]) 



2k An /C Os(t)-l 

( V 



-21n(2/iAr/) 



(80) 



Note that we have written e ^ 



: e * fcA,? and extracted the first phase 



to reconstruct the full tree order solution ^/°(r], x; k, s) 



-ikrj 



Ui(k, s)e 



ik-x 



The second identity derives from acting a d'Alembertian on (!80|) . The 
d'Alembertian passes through the tree order solution to give, 

d 2 ^°(ri,x;k,s) = ^°(r), x; k, s)d v (d v -2ik) . (81) 

Because the integrand goes like aln(a) for small a, we can pass the first 
derivative through the integral to give, 

, (\n([i 2 Ax 2 ++ ) ln(/i 2 Ax 2 _)' 



d 2 dV 



Ax 2 



Ax 2 



^°{r]',x';k,s) 



= tAn 2 ^°(r ]l x;k,s)d ri J\'y^dt^^) +21n(^)| . (82) 

We can pass the final derivative through the first integral but, for the second, 
we must carry out the integration. The result is our second key identity, 

,(]n(pL 2 Axl + ) \n([i 2 Ax 2 + _) 



d 2 dV 



V 



Ax 2 ++ Ax 2 + _ 
= i47r 2 \I> (?7,f; k,s)\ 2 In 



ty°(r]',x';k,s) 



v2a 



1+Hti)]+ [drf( 



7 i2kAr) . 



AT] 



(83) 



The final key identity is derived through the same procedures. Because 
they should be familiar by now we simply give the result, 



1 



1 



*°{rf,# ;k,s) 
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Ax 2 



iAifhr^iji, x- k, s) I dr}' e i/cA " sin(fcAry) . (84) 



/ ikArj , 



Vi 



15 



/ 
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{H 2 aa')~ l ^<9 4 
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-8 |<9 2 
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f @d 2 
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4 ^V 2 


3 


- $d 2 


6 


7 ^V 2 



Table 1: Derivative operators U^: Their common prefactor is ■ 

4 Solving the Effective Dirac Equation 

In this section we first evaluate the various nonlocal contributions using the 
three identities of the previous section. Then we evaluate the vastly simpler 
and, as it turns out, more important, local contributions. Finally, we solve 
for ^(rj, x; k, s) at late times. 

The various nonlocal contributions to ([5l) take the form, 



j ( ln(af Ax 2 



Ax 2 ++ 



ln(af Ax 2 
Ax 2 



Ax 2 



>V Arf,tf',k,s) . (85) 



The spinor differential operators U\j are listed in Table [U The constants aci 
are \i for I = 1, 2, 3, and \H for I = 4, 5. 



As an example, consider the contribution from U~y 



2 2 8 tt 4 * J \ Ax 2 ++ 



15 k 2 H 2 
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2fi 



V°{7]',x;k,s) 
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k 2 H 2 



%R-fw\r],x\k,s) x 
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3 2ii(l+i?r,) 



+ 1 



(87) 



In these reductions we have used i ^°(rj, x\ k, s) = ry°\I/ (?7, x; k, s) and 
(183|) . Recall from the Introduction that reliable predictions are only possible 
for late times, which corresponds to rj — > _ . We therefore take this limit, 



15 k 2 H 2 

2 2 8 7T 4 
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\n(n 2 Ax 2 __ 
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Coefficient of the late time contribution from each Ufj 


1 
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f {exp(2z|) + 1} 
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-z|{21n(f)-/°^( exp( - 2 ;^ ) - 1 ) 
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) 


5 


4g/°^- 2 ^'{/o" 2W ^( eXP( 7* ) " 1 ) 


+ 21n(#?/)} 
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-|z|{exp(2z|) - 1} 



Table 2: Nonlocal contributions to / rf 4 x'[S](x; x')^/°(t]', x'\ k, s) at late times. 
Multiply each term by x iHj°ty°(r), x; k, s). 



— ^iif 7 °* Q (M k,s) x ^{exp(2*-|) + l}. 

The other five nonlocal terms have very similar reductions. Each of them 
also goes to x iH r y°^°(ri, x;k, s) times a finite constant at late times. 
We summarize the results in Tableland relegate the details to an appendix. 

The next step is to evaluate the local contributions. This is a straightfor- 
ward exercise in calculus, using only the properties of the tree order solution 
(fIE|) and the fact that d^a = Ha 2 5®. The result is, 

f d V J @d 2 + — ln(oo') 9- 7 ln(oo') jpj 5 4 (x - x') ^° (?/,£'; M 

2 b 7r 2 J [ H 1 aa! 2 J 

+ y(ln(a)^f ^ln(a))^°(r/,x;fc,s) - 141n(a) ^(77, x\ k, s)j, (89) 
K " zif 7 ^°(r/, x; fc, s) x { — a - Ui^- ln(a) -2i^-\. (90) 



2 6 tt 2 ' v " ' ' ' \ 2 H w if. 

The local quantum corrections (190]) are evidently much stronger than their 
nonlocal counterparts in Table [2J Whereas the nonlocal terms approach 
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a constant, the leading local contribution grows like the inflationary scale 
factor, a = e Ht . Even factors of ln(a) are negligible by comparison. We can 
therefore write the late time limit of the one loop field equation as, 

i^^X-Xs) — ^yZiW^MM • (91) 

The only way for the left hand side to reproduce such rapid growth is for the 
time derivative in i@ to act on a factor of ln(a), 

TT 2 

iYd„ ln(a) = ^—$1 = i#a 7 ° . (92) 

We can therefore write the late time limit of the tree plus one loop mode 
functions as, 

^f^ + K 2 * 1 ^^ . (93) 

All other corrections actually fall off at late times. For example, those from 
the ln(a) terms in (l9"Uj) go like ln(a)/a. 

There is a clear physical interpretation for the sort of solution we see in 
(J9*3"j) . When the corrected field goes to the free field times a constant, that 
constant represents a field strength renormalization. When the quantum cor- 
rected field goes to the free field times a function of time that is independent 
of the form of the free field solution, it is natural to think in terms of a time 
dependent field strength renormalization, 

Tt - ? ^ ^°(i],x; k, s) , „ . , H 17k 2 if 2 , . . _ 4n 

$(ri,x;k,s) — ► — v ; where Z 2 (t ) = 1 -__l n (a) + 0(K 4 ) . 

Z 2 (t) 2e>7r 



(94) 

Of course we only have the order k 2 correction, so one does not know if 
this behavior persists at higher orders. If no higher loop correction su- 
pervenes, the field would switch from positive norm to negative norm at 
ln(a) = 2 6 7r 2 /17k 2 H 2 . In any case, it is safe to conclude that perturbation 
theory must break down near this time. 



5 Hartree Approximation 

The appearance of a time- dependent field strength renormalization is such 
a surprising result that it is worth noting we can understand it on a simple, 
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qualitative level using the Hartree, or mean-field, approximation. This tech- 
nique has proved useful in a wide variety of problems from atomic physics 
[4~3] and statistical mechanics [44] . to nuclear physics [45J and quantum field 
theory Of particular relevance to our work is the insight the Hartree 

approximation provides into the generation of photon mass by inflationary 
particle production in SQED [471 1151 139] . 

The idea is that we can approximate the dynamics of Fermi fields inter- 
acting with the graviton field operator, h^ Vl by taking the expectation value 
of the Dirac Lagrangian in the graviton vacuum. To the order we shall need 
it, the Dirac Lagrangian is [9], 



^Dirac = + -{hm^-h^ lfl id v ^-h ma ^J P ^} 



+ K 



+ K 
2 



-h 2 --h pa h 

-^^p, ct + ^"p^^ + ^(^^p)^ + ^">ma-]^V^^ + 0(k 3 ). (95) 



Of course the expectation value of a single graviton field is zero, but the 
expectation value of the product of two fields is the graviton propagator 
[131, 



{n\T[h txv {x)h pa {x')\ I n) 

= iA^(x; x' 



rpA 

V-v 1 pa 



+ iA B (x; x' 

The various tensor structures are, 

2 



rpA 

pv 1 pu 



2r]p(pVa)v 



D-3 



VpvVpa 



rpB 

t^ /1 pa 



rpB 

V-v 1 pa 



+ iA c (x: x') 



rpC 

Vv 1 pa 



rpC 

fw 1 pa 



(D-2)(D-3) 



(D-3)M + ^ (D-3)6X + ri P a 



0x0 



(96) 

(97) 
(98) 



Parenthesized indices are symmetrized and a bar over a common tensor such 
as the Kronecker delta function denotes that its temporal components have 
been nulled, 

r v = %-8ffi , ^ = ^ + ^°. (99) 

The three scalar propagators that appear in (I9"6"j) have complicated ex- 
pressions which we omit in favor of simply giving their coincidence limits and 
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the coincidence limits of their first derivatives |5Uj . 



lim iAa(x; x') 

X 1 — >x 

lim dniA.^(x; x') 

x'— *x 

lim iAsix; x') 

x'— >x 

lim dniA B (x; x') 

x'— >x 

lim iAc(x; x') 

x'^x 

lim d^iAcix; x') 



H D - 2 Y{D - I) 

(47r)f r(f) 

H D-2 _ Vj 

(47r)f r(f) 
H p-2 r ( D _ Vj 

(47r)f r(f) 

= lim d'iAnix: x') 



7T 



-vrcot^-DJ +21n(a) \, (100) 



x Ha5l = lim d'iA A (x;x') , (101) 



x 



D-2 



x 



^-2 r (£) _ 1) 

(4tt)t T(f ) ~ (D-2)(D-3) 
= lim 9^Ac(x; x') . 



(102) 
(103) 

(104) 
(105) 



We are interested in terms which grow at late times. Because the 5-type and 
C-type propagators go to constants, and their derivatives vanish, they can 
be neglected. The same is true for the divergent constant in the coincidence 
limit of the A-type propagator. In the full theory it would be absorbed into 
a constant counterterm. Because the remaining, time dependent terms are 
finite, we may as well take D = 4. Our Hartree approximation therefore 
amounts to making the following replacements in (I95p . 



III 

471-2 

H 2 



ln 0) [w»+V^r 2 VV 

VppVua'^Vp.aVup^^'VpuVpiT 



(106) 
(107) 



It is now just a matter of contracting (I106H107P appropriately to produce 
each of the quadratic terms in (1^5]) . For example, the first term gives, 



Kj _ r>— — 



K 2 H 2 

2%2~ 
k 2 H 2 

2 5 7T 2 



ln(a) 
ln(a) 



3 + 3 



rinpHua + VpvVvp - ZVpvVpv ^ifiV, (108) 

. (109) 



The second quadratic term gives a proportional result, 

2 _ K 2 H 2 
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-h pa h p<j m^ 



2% 2 



ln(a)[9 + 3-6 mftH . 



110) 
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The total for these first two terms is ~^? 2 ln(a)^fi <ft^ . 

The third and fourth of the quadratic terms in (1951) result in only spatial 
derivatives, 



-k 2 H 2 

~2^ 2 
3k 2 H 2 



8 



ln(a 
ln(a) 



1 + 1-6 
3+1-2 



, (111) 
v^jM . (112) 



The total for this type of contribution is ln(a)^i^^. 

The final four quadratic terms in (195]) involve derivatives acting on at 
least one of the two graviton fields, 



-k 2 H 2 
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2H 2 
k 2 H 2 
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K 



2 6 7T 2 
K 2 H 2 

2%2" 
. 



Ha 
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3 + 1-2 



fj^-fj^ , (H3) 
^7°^*, (114) 
3 + 1-2Jt7 w *7^ *, (115) 

(116) 



The second of these contributions vanishes owing to the antisymmetry of the 



Lorentz representation matrices, J^ u = |[7 M ,7^], whereas rj „ p 7 Al </ p0 = — %7° 



Hence the sum of all four terms is 
Combining these results gives, 
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167T 2 



ln(a)^^+0(K 4 ).(118) 



If we express the equations associated with (11181) according to the perturba- 
tive scheme of Section 1, the first order equation is, 



k 2 H 2 - k 

i^K A ^ L {r],x]k,s) = 1 ^^iH^°{r],x-,k,s)[24a-28i—ln(a)j . (119) 



2 6 7T 2 
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This is similar, but not identical to, what we got in expression (1901) from the 
delta function terms of the actual one loop self-energy fl3]). In particular, the 
exact calculation gives ya — 14i-4 m(a), rather than the Hartree approxima- 
tion of 24a — 28i-4 ln(a). Of course the ln(a) terms make corrections to 
which fall like ln(a)/a, so the real disagreement between the two methods is 
limited to the differing factors of y versus 24. 

We are pleased that such a simple technique comes so close to recover- 
ing the result of a long and tedious calculation. The slight discrepancy is 
no doubt due to terms in the Dirac Lagrangian (1951) which are linear in the 
graviton field operator. As described in relation (132|) of section 2, the lin- 
earized effective field ^i(x; k, s) represents a~~^~ times the expectation value 
of the anti-commutator of the Heisenberg field operator ipi{ x ) with the free 
fermion creation operator b(k,s). At the order we are working, quantum 
corrections to ^(x; k, s) derive from perturbative corrections to ipi{x) which 
are quadratic in the free graviton creation and annihilation operators. Some 
of these corrections come from a single hhipip vertex, while others derive 
from two hipip vertices. The Hartree approximation recovers corrections of 
the first kind, but not the second, which is why we believe it fails to agree 
with the exact result. Yukawa theory presents a fully worked-out example 
0, E, El] in which the entire lowest-order correction to the fermion mode 
functions derives from the product of two such linear terms, so the Hartree 
approximation fails completely in that case. 

6 Discussion 

We have used the Schwinger-Keldysh formalism to include one loop, quantum 
gravitational corrections to the Dirac equation, in the simplest local Lorentz 
and general coordinate gauge, in the locally de Sitter background which is 
a paradigm for inflation. Because Dirac + Einstein is not perturbatively 
renormalizable, it makes no sense to solve this equation generally. However, 
the equation should give reliable predictions at late times when the arbitrary 
finite parts of the BPHZ counterterms are insignificant compared to the 
completely determined factors of ln(aa') on terms ([7}l9]) which otherwise have 
the same structure. In this late time limit we find that the one loop corrected, 
spatial plane wave mode functions behave as if the tree order mode functions 
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were simply subject to a time-dependent field strength renormalization, 

Z 2 (t) = 1 - — GH 2 \n(a) + 0(G 2 ) where G = IQtxk 2 . (120) 

If unchecked by higher loop effects, this would vanish at ln(a) ~ 1/GH 2 . 
What actually happens depends upon higher order corrections, but there is 
no way to avoid perturbation theory breaking down at this time, at least in 
this gauge. 

Might this result be a gauge artifact? One reaches different gauges by 
making field dependent transformations of the Heisenberg operators. We 
have worked out the change (1681) this induces in the linearized effective field, 
but the result is not simple. Although the linearized effective field obviously 
changes when different gauge conditions are employed to compute it, we 
believe (but have not proven) that the late time factors of In (a) do not 
change. 

It is important to realize that the 1PI functions of a gauge theory in a 
fixed gauge are not devoid of physical content by virtue of depending upon 
the gauge. In fact, they encapsulate the physics of a quantum gauge field 
every bit as completely as they do when no gauge symmetry is present. One 
extracts this physics by forming the 1PI functions into gauge independent 
and physically meaningful combinations. The S-matrix accomplishes this in 
flat space quantum field theory. Unfortunately, the S-matrix fails to exist 
for Dirac + Einstein in de Sitter background, nor would it correspond to an 
experiment that could be performed if it did exist [521 [531 [51] . 

If it is conceded that we know what it means to release the universe in a 
free state then it would be simple enough — albeit tedious — to construct an 
analogue of ipi(x) which is invariant under gauge transformations that do not 
affect the initial value surface. For example, one might extend to fermions 
the treatment given for pure gravity by [55] : 

• Propagate an operator-valued geodesic a fixed invariant time from the 
initial value surface; 

• Use the spin connection A^ cd J cd to parallel transport along the geodes- 
ic; and 

• Evaluate ip at the operator-valued geodesic, in the Lorentz frame which 
is transported from the initial value surface. 
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This would make an invariant, as would any number of other constructions 
[56] . For that matter, the gauge-fixed 1PI functions also correspond to the 
expectation values of invariant operators [57]. Mere invariance does not guar- 
antee physical significance, nor does gauge dependence preclude it. 

What is needed is for the community to agree upon a relatively simple 
set of operators which stand for experiments that could be performed in de 
Sitter space. There is every reason to expect a successful outcome because 
the last few years have witnessed a resolution of the similar issue of how to 
measure quantum gravitational back-reaction during inflation, driven either 
by a scalar inflaton [581 EE EE EE] or by a bare cosmological constant [62J . 
That process has begun for quantum field theory in de Sitter space [53" 1 15"4" l 156] 
and one must wait for it to run its course. In the meantime, it is safest to 
stick with what we have actually shown: perturbation theory must break 
down for Dirac + Einstein in the simplest gauge. 

This is a surprising result but we were able to understand it qualitatively 
using the Hartree approximation in which one takes the expectation value of 
the Dirac Lagrangian in the graviton vacuum. The physical interpretation 
seems to be that fermions propagate through an effective geometry whose 
ever-increasing deviation from de Sitter is controlled by inflationary graviton 
production. At one loop order the fermions are passive spectators to this 
effective geometry. 

It is significant that inflationary graviton production enhances fermion 
mode functions by a factor of ln(a) at one loop. Similar factors of ln(a) have 
been found in the graviton vacuum energy [63, 64J. These infrared logarithms 
also occur in the vacuum energy of a massless, minimally coupled scalar with 
a quartic self- interaction [651 EE], and in the VEV's of almost all operators 
in Yukawa theory [5T] and SQED [67J. A recent all orders analysis was not 
even able to exclude the possibility that they might contaminate the power 
spectrum of primordial density fluctuations [68] ! 

The fact that infrared logarithms grow without bound raises the exciting 
possibility that quantum gravitational corrections may be significant dur- 
ing inflation, in spite of the minuscule coupling constant of GH 2 < 1CT 12 . 
However, the only thing one can legitimately conclude from the perturbative 
analysis is that infrared logarithms cause perturbation theory to break down, 
in our gauge, if inflation lasts long enough. Inferring what happens after this 
breakdown requires a nonperturbative technique. 

Starobinskii has long advocated that a simple stochastic formulation of 
scalar potential models serves to reproduce the leading infrared logarithms of 
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these models at each order in perturbation theory [69|. This fact has recently 
been proved to all orders [701 EQ- When the scalar potential is bounded 
below it is even possible to sum the series of leading infrared logarithms and 
infer their net effect at asymptotically late times [72]! Applying Starobinskh's 
technique to more complicated theories which also show infrared logarithms is 
a formidable problem, but solutions have recently been obtained for Yukawa 
theory [M] and for SQED [57]. It would be very interesting to see what 
this technique gives for the infrared logarithms we have exhibited, to lowest 
order, in Dirac + Einstein. And it should be noted that even the potentially 
complicated, invariant operators which might be required to settle the gauge 
issue would be straightforward to compute in such a stochastic formulation. 



7 Appendix: nonlocal terms from section 4 



It is important to establish that the nonlocal terms make no significant con- 
tribution at late times, so we will derive the results summarized in Table HJ 
For simplicity we denote as [U 1 ] the contribution from each operator in 

Table[TJ We also abbreviate \1/ (77, x; k, s) as ty°(x). 

Owing to the factor of I /a' in U^, and to the larger number of derivatives, 
the reduction of [U 1 ] is atypical, 



[IT 1 ] 



K 1 



2 8 tt 4 a ' 



/ <Fx'- 



, 1 [ ln(/x 2 A;c 2 + ) ln(/z 2 Aa^ J 



V 



>y°(x^ 



Ax 2 ^ 



Ax 2 



;i2i) 



-IK 



2 6 vr 2 a 



-IK 



2 6 vr 2 a 



y°y°(-2ik + d 7 



{rri p 2ikAri _ i 

+«9 2 rdr]'{-2Hn')ln(2fxAr]) 



122) 



(r) 



-) 2 



(2ik - H)e 2ik ^ + ^ 3H 2 



k 2 H 2 
2H 2 



V + JI 
2 



+ 



2H 3 r] 



1 + Hrj) (1 + HriY 



(123) 



. _|_ / / 



[1 + H V ) 3 

5 - Aikrj - ^ 



[1 + Ht]) 2 



+ 



6? A 
H 



[1 + Ht]) 2 1 + Hi] 



■ (124) 



25 



This expression actually vanishes in the late time limit of rj — > . 
[U 2 ] was reduced in Section 4 so we continue with [U 3 ], 
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Ufj has the same derivative structure as Ufj, so \U A ] follows from (11281) . 

k 2 # 2 n ,„ /4 L r,- „ ,i n . ,/e 2ifcA "-lxl 



^8fc7 $°(x)j21n[(l + #77)] + f drf( 



^piH-fVlz) x 8,| / V( 
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C/^- has a Laplacian rather than a d'Alembertian so we use identity flHUl) 
for [£7 5 ]. We also employ the abbreviation kAr] = a, 

m = - W> , (132) 
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, , f 2a , sin(t) . , , r f 2a , /cos(t)-l\ „, 
x^-cos(a)/ dt — ^ + sin(a) / dt( ^ +21n(— M133) 
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2 b 7T 2 ii Jrii l J0 V £ 7 J 

£/f„- has the same derivative structure as Uf, but it acts on a different inte- 



26 



grand. We therefore apply identity (I8"4"j) for [U 6 ], 

I^^V^v{^-^-}*V), (136) 

2 7 7 2 

= 7^—$V 2 x (-^Tr 2 )*- 1 * ^) [ dr]'e ikAl, sin(kAri) , (137) 
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